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Abstract 

In this paper the equation of state of hot neutron matter is calculated. 
Involving the Oppenheimer-Volkov- Tolman equation global parameters 
of a neutron star at the finite temperature are obtained. The objective 
of our work was to study the influence of the temperature on the main 
parameters of a neutron star. 



1 Introduction 

This paper is concerned with a neutron star and its macroscopic parameters 
such as the mass M and the radius R which are influenced by the temperature. 
Both the mass and the radius as well as the cooling evolution are determined 
first of all by the equation of state. 

Considering the matter of a neutron star one should take into account not only 
neutrons but as the most elementary model has it neutrons, protons and leptons. 
This paper presents a basic model of neutron star matter including interactions 
among nucleons in the Hartree approximation Increasing interest in neu- 

tron matter at finite temperature has been observed recently in relation to the 
problems of hot neutron stars and of protoneutron stars and their evolutions in 
particular. Theories concerning protoneutron stars are being discussed in works 
by Prakash et. al. ||. Some other elements like hyperons, mesons or quarks 
could be also found in the interior of a neutron star but their relevance is not 
going to be included in our work. 

One can divide this paper into two parts. In the first one thermodynamic proper- 
ties of Fermi gas, which consists of neutrons, protons and leptons, are examined. 
The properties of the quark matter at finite temperature were also looked into 
in [Q. The range of temperatures considered vary from 0-50 MeV. Analytical 
forms for the pressure, energy density and fermion number density has been 
calculated in order to solve the OTV equation. Solution of this equation is the 
main subject of the second part. 
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2 General Theory 

The aim of this paper is to present a rudimentary approach to the equation of 
state for neutron stars at the finite temperature. In such an approach the neu- 
tron star matter consists of electrically neutral plasma which comprises protons, 
neutrons and electrons. The Lagrange density function in this model is given 
by 



C = ii/jj^d^tp — ipMip + v< ipi/j > ipip — -v < i/iip > 2 + 



(1) 



1 



2 2 

* Yl L?f d n L t - 9f(LfHe Rf + h.c.) - j-R 
/=i /=i 



where k — 87rG/c 4 and R is the Ricci curvature scalar. The fermion fields are 
composed of neutrons, protons and electrons, muons and neutrinos 



, L 2 = 
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and the Higgs field H takes the form of 

1 



H = 



V2\V 



(3) 



Nucleon masses are given by 



M 



m n 

TO 



V 



(4) 



The electron and muon masses equal m e — giV/V2, — g-2Vj\f2 respectively 
where V = 240 GeV. The model describes nuclear interaction in the Hartree 
approximation for » ^ Q. Introducing the interaction one can obtain the 
effective nucleon mass 

M* = M — v < Tpip > (5) 



The constant term appearing in (|l|)is responsible for the negative pressure 

1-9 

P = -v < tpip > z 



(6) 



This model is the simple approximation of the relativistic mean filed theory |T| . 
The equation (||) is a highly nonlinear equation due to the average < iptp > 
which is defined as 
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Figure 1: The effective nuclear mass 

The solution of the equation (||) help us to define the parameter S = m e ///m. 
The term m e // denotes fermion effective mass which for nucleons is written as 
m e ff = M* whereas for leptons m e ff = m thus 6—1. Fig. 1 shows the relation 
between the parameter 5 and the dimensionless Fermi momentum x for three 
different temperatures T = 0, 10, 50 MeV. The simplest case corresponds to 
L2 parameter set |Q. In the relativistic mean filed theory it corresponds to 
m s — 520 MeV and g s = 10.47 what gives 



The fermion number density, energy density and pressure inside the star are 
defined locally as quantum averages. One can calculate such quantum averages 
employing the following equation 



< A>=±Tr(e- fj(n -^f^ N ' ) 



Z 



A) 



(7) 
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where A is an observable, H is the Hamiltonian of the system and 

Z = Tr(e- f3{n -^f^ Nf) ) 

is a partition function. In this equation /3 = 1/fesT. Let us now define some 
operators indispensable in our calculations: the fermion number operator 

N f = EK„, f c Kaj - 4,./w) (8) 

where the index / stands for neutrons, protons and electrons (/ = p,n,e) and 
c k / ' c k a f anc ^ f' a f are crea tion and annihilation operators for par- 
ticles and antiparticles respectively, the fermionic Hamiltonian Ji (ignoring the 
zero-point energy) 

n XX'k.,;<- k .,../-k,., + d L,A.,/) (9) 

/ k . ct 



with ejj. j = c^Jfi 2 k 2 + rrijC 2 and the pressure with an isotropic distribution of 
momenta given by 

P f = I E E l>'^'- k ,,r'-k,,r + (10) 
/ k,<r 

where velocity equals «k / = ^c 2 / e k /• ^he mean number of fermions is de- 
termined by the following equation 

< c Lj c k,«j >= exp(/3( £kj - M/ )) + l (U) 

where jif stands for the fermion chemical potential. Neutrons, protons and 
electrons are in /3-equilibrium which can be described as a relation among their 
chemical potentials 

flp + H e = fl n (12) 

where fi p , /i„ and fx e stand for proton, neutron and electron chemical potentials 
respectively. If the electron Fermi energy is high enough (greater then the 
muon mass) in the neutron star matter muons start to appear as a result of the 
following reaction 

e~ — > n~ + v e + T7JI (13) 

The chemical equilibrium between muons and electrons can be described by the 
condition 

/V = Me (14) 



The equation (12) together with the charge neutrality n e + — n p allows us 
to determine the equation of state in terms of only one parameter x n . It defines 
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the dimensionless Fermi momentum of a neutron x n = (Tik n ) / (m„c) . Using the 
stated above conditions and the equation 

M / = m / c 2 (l + 4) 1 / 2 (15) 

the relations between x n = x and x p can be obtained 



In the macroscopic limit an integral is allowed to replace a sum 

V 



£ - tLs I d3k (16) 



k 

and for each fermion the quantum averages n f representing the particle number 
density n/ =< Nf > /V, £f the energy density and Pf the pressure can now 
be written as 

n / (M/ ) T)=n e(r,y) (17) 

effaf, T) = c 2 p( H ,T) = e x(r, y) (18) 

P(fx f ,T) = P ^(r,y) (19) 

The forms of these quantum averages are determined by the functions Q(r,y), 
x( r ,y) and <J>(r, y) which are presented below. 



z^dz 1 



a/z 2 + S 2 y 2 exp {yj z 2 + S 2 y 2 — ry) + 1 

-} 



exp (y/ z 2 + S 2 y 2 + ry) + 1 ' 



exp (y z + o z y z — ry) + 1 



( 



exp ( y/ z 2 + S 2 y 2 + ry) + 1 



1 f°° 1 

x (r,y) = — / dzy/z 2 +5 2 y 2 { (22) 

T Jq exp (y/z 2 + S 2 y 2 - ry) + 1 

+ 1 , 

exp {y/z 2 ^^ 2 ^ 2 + ry) + 1 

with the following dimensionless variables S = m e ff/m, y — mc 2 /ksT, z — 
Tikc/hsT and r = /if /mc 2 . The forms of the dimensionless variables are the 
same as those used by Weldon in his work ||, m. The fermion chemical potential 
Hf is given by the relation (|l5|). Mutual relations between y and temperature 
as well as between r and /i/ result in dependence of functions 0, $ and \ on 
the dimensionless Fermi momentum x and the temperature T. 
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When the temperature equals zero the Shapiro for result free nucleons can be 
reproduced || 

°) = 7 ^t{xVTTx 1 {2x 2/3 - 1) + \n{x + y/l + x*)} (23) 

which in the nonrelativistic limit x << 1 yields 

SM)-^* 6 (24) 

and the equation of state takes the form of the polytrope T — 5/3. The forms 
of the functions 0, \ an d indicate their relations with the functions H n (r, y) 
and G n (r,y) which are used in order to evaluate thermodynamic properties of 
the matter 



H„(r,y) = / —===={ (25) 

( 



exp(-\/ z 2 + (5 2 y 2 + ry) + 1 



G„(r, y) 



Jo 



exp(-\/ z 2 + (5 2 y 2 + ry) + 1 



exp(y / z 2 + 5 2 y 2 — ry) + 1 



} 



(26) 



The case where the interaction in the Lagrange function ( |l[) i s neglected (free 
nucleons) is equivalent to the one where S in equations ( [25| , p6"l ). These two 
terms in both equations ( 25 2^) correspond to the contribution of particles and 
antiparticles, respectively and the functions H n (r, y) and G n (r, y) can be written 

as 



H„(r,y) = h n (r,y) +h n (-r,y) 
G n (r,y) = g n (r,y) - g n (-r,y) 



(27) 



Both the pressure P/ and the energy density e/ of the fermion system can be 
expressed with the use of the functions presented above 

Pf = ^ Po ^ 4{H5 ^ y P^ + H 5( r n,yn) + H 5 (r e ,y e )} (28) 

where Pq = (m/c 2 )/A 3 , the fermion Compton wavelength A = h/(rrifc) and 
A T = (ch)/(k B T) 

£f = -^P (-^f{H 5 (r p ,y p )+H 5 (r n ,y n ) + H 5 (r e ,y e ) (29) 

IT AT 

+ H 3 (r n ,y n ) + H 3 (r p ,y p ) + H 3 (r e ,y e )} 
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Our first step is to find the pressure and energy density for nucleons and muons 
(6=1) which corresponds to the nonrelativistic limit of the functions G n (r, y) 
and H n (r, y) . Such a limit means either the case of large mass or low tempera- 
ture. Introducing the new variable u> = exp (y — x 2 + y 2 ) the following forms 
of the functions are achieved 

*(r, y) = I [ - \ y \ s l )+(J (30) 

, . f 1 (-\nuj) 3 / 2 (2yS-lnuj) 3 / 2 {y5-lnuj)duj , 
9dr, V) = J o e vV-r )+u ( 31 ) 

The calculation of these integrals has been performed on the basis of the method 
presented by Weldon in his work H, |(|. Making use of the fact that | In ^\ < 1 
in this approximation and expanding the numerators under that assumption, the 
functions h§ and g$ can be written in the form of 

Mr, V) = ^ E {2dy)kT{ l_ m I eVls l )+u (32) 

( \ (9 na/»f f 1 (-lnc) 3 / 2 (2^-ln^/^ 

Having integrated the obtained equations term by term the final forms of the 
function emerge: 

h 5 (r,y) = (2 y S)^T( 5 -)± I^^^Li^-^) (34) 



S 1 1 

S .Ar,y) = ( art V^_ ) g_ f _ ( _). ( 35) 

(tfU M M-^"-* > )nl + k > + rl -l + l>)Li mk (-e"< r -»)) 

In order to calculate the contribution of electrons to the total pressure Pf and 
the energy density £/, the relativistic case should be considered. In this very 
case the variable y tends to zero. It is necessary to calculate the mentioned 
above functions H n (r,y) and G n (r,y) in the relativistic limit, where again the 
method used by Weldon in Q, J(| has been involved. These functions obey the 
recursion relations 



dGn+l V „ , V 2r tt 

H ti n -i 
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(36) 
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Thus the functions Hi (r, y) and G\ (r, y) together with the initial conditions 
G„(0,0) = 

ff„(0,0) = 2(l-2 2 -")r(n-l)C(n-l) (37) 

are sufficient to determine the functions H^(r,y) and Hs(r,y) which in turn are 
indispensable to calculate the pressure and energy density. The identity 

1 1 00 y 

expy + 1 = 2 ~ 2 ^y 2 + 7r 2 (2n+l) 2 ^ 

originated from Dolan and Jackiw || gives us the possibility to write the func- 
tions Hi and G\ as 

r°° dz 
Jo ^{z 2 + y 2 ) 

[z 2 + y 2 (\-r 2 ) + {2n+l) 2 ir 2 ]dz 



Gi(r,y) 



oo „ 

-4ry^ / 



[z 2 + y 2 (l - r 2 ) + (2n + 1) 2 tt 2 ] 2 + Ait 2 r 2 y 2 {2n + l) 2 
[z 2 + y 2 (\-r 2 )- (2n+l) 2 Ti 2 ]dz 



[z 2 + y 2 (l - r 2 ) + (2n + 1) 2 tt 2 ] 2 + 47r 2 r 2 y 2 (2n + l) 2 

(40) 

The integrands in these equations are multiplied by the convergent factor z e 
and after its expansion as a power series in y the term by term integration is 
performed. In the next step the summation over n is carried out. In the last 
stage the limit e — ► enables us to obtain the final result which in the first 
approximation achieves the form of 

ffi(r,y) = -( 7 + ln(-)) (41) 

G 1 (r,y)=ry (42) 

Knowing the form of functions Hi (r, y) and G\ (r, y) and the initial conditions 
( |37| ) it is possible to calculate the functions H^{r,y), H$(r,y), Gz{r,y) 

H z {r lV ) = l -{2r 2 + 2 7 - l)y 2 + -^lnfoAO (43) 

H 5 (r,y) = ^(8r 4 - 24r 2 - 12 7 + 9)y 4 - ^y 4 lny/n (44) 

G 3 (r,y) = ^r(2r 2 -3)y 3 (45) 

These functions are employed to express the electron pressure and energy density 
according to the equations (28) and (|30|). Consequently the relation between 
the total pressure Pf being the sum of the leptons and nucleons pressures and 
the dimensionless Fermi momentum x is presented in Fig. 2. The curves in 
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Figure 3: The energy density x{ x ) as the function of the Fermi momentum 

x = pp/mc 



Fig, 2 are parameterized by the temperature and 8. The curves for the zero 
temperature limit are obtained for two cases with and without (v = 0) the 
presence of interaction. The curve for the temperature T = 50 MeV is the case 
with interactions. 

Analogous description can be made analyzing the connection between the 
energy density \ an d the parameter x. See Fig. 3. 



3 The neutron star 

The most important factor determining the structure of a neutron star is the 
equation of state. This very equation makes possible to describe a static spher- 
ical star solving the OTV equation. 

^1 - 47rr 2 p(r) (46) 



10 



dr r 2 C 2 Cl 2Gm(r) x v ' 

Having solved this equation the pressure -P(r), mass m(r) and density p(r) were 
obtained. To achieve the total radius R of the star the fulfillment of the condition 
P{R) = is necessary which allows to determine the total gravitational mass 
of the star M(R). 

Introduction of the dimensionless variable £ which is connected with the variable 
r by the relation r = a£ (a = 1 km) enables us to define the functions P(r), 
p{r) and m(r) in the following form 

p(r)=PcfHZ) (48) 

P(r) = P c u(0 (49) 
m(r) = M u(0 (50) 

Some more parameters, namely 

^ = ^f- 52 

and 

are also needed to achieve the useful form of the OTV equation 

% = -W0*+™®) V ® +bmt{ S? (54) 



dv 



with 



rfiOH 2 (55) 

r a = -^r 1 (56) 



These equations can be solved specifying the central neutron energy density 
being the energy density for r = 0. 

The equation of state is the function of the temperature and the neutron chem- 
ical potential which changes with the radius r. Therefore the changes of the 
radius influence the parameters of the star. Using the results obtained in pre- 
vious chapter it is now possible to write the functions P(r) and p(r) in the 
form 

p(n,T) = pox{x,T)=p c fi{Q (57) 
P $(x,T) = P c u(0 (58) 
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Figure 4: The R-M diagram for the neutron star with interaction and without 
(v = 0). 

and the variable x can be obtained 

x = X - 1 (-fH0) (59) 
Po 

Pc u(0 = PoHx- 1 (-f i ),T) (60) 
Po 

The function u(£) takes the form 

«(0 = §*(x" 1 (^/ § (0).r)- (6i) 

The solution of the Oppenheimer-Volkoff-Tolman equation depicts the mass- 
radius relation. Fig. 4 allows us to compare the zero-temperature mass versus 
radius relation with the other temperature cases with and without (v = 0) 
interaction . 

Fig. 5. shows the changes of the radius as a function of the neutron density 
p c in the center of the star. 
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Figure 5: The neutron star radius R(km) as a function of the neutron density 
p c in the center of the star. 
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4 Conclusions 



Herein presented solutions of the Oppenheimer-Volkoff-Tolman equations of hy- 
drostatic equilibrium concern hot neutron stars. The neutron star matter in 
this model consists of neutrons, protons and leptons (electrons, muons) being 
in j3- equilibrium under the assumption that the temperature is different from 
zero. It varies from to 50 MeV. The objective of our work was to study the 
influence of the temperature on the main parameters of a neutron star. In order 
to achieve the proper form of the equation of state, which is determined only 
by the neutron Fermi momentum, it is necessary to calculate either the low 
temperature or the high temperature expansion of the integrals (??) and (g7j). 
This very simple model presents a few global properties of the hot neutron star 
such as the mass and the size. The parameters of neutron stars obtained in this 
simplified model considering zero temperature and cases of finite temperatures 
vary from one another. In each of the mentioned above cases stars with the 
same energy density inside are considered. However, their baryon numbers are 
different which makes each of them a different star with specific baryon num- 
ber. In this situation we do not deal with the thermal evolution of one star 
with conserved baryon number but several separate cases. The star whose pa- 
rameters at the temperature of 50 MeV are as follows: for 8 = 1 M = 0.62 M Q 
R = 8.08 km, for 8 ^ 1 M = 1.3 M©, R = 8.00 km At zero temperature limit 
the star is characterized by the mass 0.7 Mq and the radius 10.54 km for 8 = 1 
and the mass 1.24 M© and radius 7.98 km for 8^1. It is obvious that some 
more extended models e.g. those including boson fields should be examined. 
Therefore we would like to continue the subject in our next papers 
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